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We consider the Cauchy problem of the gravitational wave with the initial data 
distributed only around the brane in the one brane model of Randall and Sundrum, 
and examine its behavior as t oo . Then we find its leading behavior is unlike an 
ordinary flat 5-dimensional space-time. Such a signal shows that the Huygens principle 
is violated on the 4-dimensional brane world even asymptotically and also shows the 
difference between compact and non-compact brane worlds. Some comments are also 
given related to AdS/CFT correspondence . 
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§1. Introduction 



Recently, inspired by the string theory,0) there are many works called as brane world 
scenario,&'&'&'B which claims "our world is confined to the 4-dimensional brane that is 
embedded in higher dimensional space time" . With respect to the existence of higher di- 
mensional space-time, this idea resembles the old Kaluza-Klein compactification scenario, 
but has great difference from it that the wave function does not spread over the extra di- 
mension but has its support only on the brane. Various papers utilize such a 5-distributional 
wave function to explain the hierarchy or others. Around such a work, one of the most 
surprising claim is that we do not need the compactification to get the 4-dimensional gravity 
near the brane; this is "An alternative to compactification" by Randall and Sundrum^ which 
we refer to as RS2 model henceforce. So if you accept the extra-dimension, then there are 
three ways to utilize it, an old Kaluza-Klein, brane worlds with compact extra dimensions 
and with non-compact ones. But, then, "What is the difference between them?", or "Cannot 
we see the difference except at the high energy?" To answer these questions, we examine the 
Cauchy problem of the gravitational wave with the initial data distributed only around the 
brane in the RS2 model, that is the third possibility of the extra dimensions. Our restriction 
to the initial data is very natural in the brane world scenario since we are confined on the 
brane and any physical phenomena takes place locally, so if we disturb the space-time then 
such a signal is distributed only around the brane. 

The organization of this paper is as follows. In the next section, we review the RS2 model 
and formulate the problem which we want to solve. In the third section, we solve and analyze 
the above Cauchy problem. Our main result is the tail from the (5-distributional source 
of gravitational wave. RS2 model is also interesting from the fact that it has holographic 
dual interpretation. We also make comments on the above behavior from the 

holographic viewpoint. In the final section, we give some comments and discussions. 

§2. Setting the problem 

As discussed in the introduction, we want to solve the Cauchy problem of the gravitational 
wave in RS2 model. So we first review RS2 and set the problem appropriately. We will mainly 
follow the framework given by Ref. IT2|). 



In the RS2 model, we start with the five-dimensional manifold (and its mirror M^) 
with two-boundaries, one is the brane boundary Ei which is introduced by hand and another 
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is the infinite far-away boundary E2 from the brane. So the Lagrangian is given by 
1 
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[2-1) 



J J El JdM^ 

We write the Lagrangian on the half of RS2 geometry. Actually an orientation reversed 
five-dimensional manifold must be glued at the brane Ei. In k,^ is the unit 

length of 5-dimensional bulk, is a 5-dimensional cosmological constant, g is the induced 
metric on the boundary, the second term is the matter Lagrangian and the third term 
represents the extrinsic curvature contribution from the boundary dM^ = i7i U called a 
Gibbons-Hawking term.© The necessity of the third term is due to the requirement that, 
in the variational principle we only demand that the infinitesimal variation itself vanishes 
uniformly along the infinite far-away boundary, but the derivative of it does not vanish on 
the boundary.0) Under such a requirement the extrinsic curvature term is canceled with 
the contribution of the derivative of variation from the bulk. However in the Lagrangian 
( |2-1| ), we have also the brane boundary term and we want non-zero fluctuation on the brane, 
so above statement does not hold on the brane side. Under the appropriate infinitesimal 
variation, the Lagrangian (|2-1| ) changes as. 



JMs ^ ^ JT,i 

-[ d''x^g5g^''{K^,~~g^,K), (2-2) 



SO that the first term leads to the equation of motion and the second and third terms leads 
to the Israel junction conditionlll* due to the fact 5g^'^ 7^ on the brane: 



G^u — -g^uA^ — 0, 
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AK^, = T--"- - -g,.T^^''- . (2-3) 

Here we use AKfj_i, as the difference of extrinsic curvature across the brane. The apearance of 
the differnce AKfj_i, is due to the another contribution of the third term in ( p-2[ ) from mirror 
manifold M^. This Israel junction condition with respect to K/j^^, means that the energy and 
the pressure on the brane is balanced against the curvature of the bulk. So when the brane 
term consists only of the tension term £matter = A, we are able to get the RS2 geometry 
which preserve the 4-dimensional Lorentz symmetry, 

ds^ = e^-^'''~^\'-''\^/'dxl + dr\ 
e~^"°^^ A , 12 
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if we require the Z2-orbifold condition K^i,{r = ro + e) = —K^^ij = ro — e) . Here we select 
the coordinate patch under which the brane is located at r = tq. 

Then we consider the fluctuation around the above geometry. If we vary the metric 
ds^ ds^ + hfj^^dx^dx" and require that this is also the solution of the equation of motion 



with the boundary condition ( p-3|) , then we get the equation of motion for the gravitational 
wave, 

V^'dMpX + ^ [dl - jd,^ h,x = 0, 

dzhpx = aX z = d. (2-5) 



Here we consider the mode satisfying h^^ = d^h^^ = h^^ = h^^ = and have used the change 
of variables exp (— 2r/Z) = 1/5^, exp (— 2ro//) = h^^^, = h^y/z^ in the region r > tq. In 
this expression, we easily see that the mode hp\{x^ , x^, x^, x^) independent of 5-th coordinate 
5, obeys the 4-dimensional massless equation of motion, and in terms of the original h^^, 
such a mode is of the form h^^j oc exp (— 2r/Z) so that the gravity is "localized" around the 
brane. However, at this stage, we wonder why such a definite form of gravitational wave can 
be created by the four dimensional localized people. In generic situation, when we create 
a disturbance just on the brane, it is physically natural to think that the created mode 
is not such a harmonic mode on the 5-dimensional space-time but a localized disturbance 
around the 4-dimensional braneQ. So if we consider the interaction with the gravitation, 
we must consider not only the above harmonic stationary wave, but also the general initial 
condition distributed around the brane and solve its asymptotic behavior to compare with 
the 4-dimensional pure gravity. So our problem is to solve ( |2-5D under the following initial 
condition, 

h^u(t 0, X, z) ffiui^'^j ^) y 
dth^yit = 0, X, z) = kf,u{x, z). (2-6) 

where / and k has its support around the brane. Of course, in principle, we must consider 
how the gravitational wave is created by the localized matter on the brane and then solve 
the above equation, but in this paper we don't consider such a problem for simplicity. 

*' This is like the situation under which we hve on the earth and make a disturbance on the ground in 
various ways. Such a wave form of disturbance depends on how we create it, the colhsion of a meteorite 
or the walking. In comparison with such a situation, the above context with localized gravity looks like a 
Rayleigh wave's one a little. 
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§3. Analysis of asymptotic behavior 

In the previous section, we have formulated the Cauchy problem of gravitational wave in 
RS2 model. In this section, we examine the asymptotic behavior of the solution as t ^ oo. 
For this purpose, we recall the way how to solve the Cauchy problem in general. When we 
consider the Cauchy problem of some self-adjoint positive definite operator A, 

—<P{x,t) + A^{x,t)=0, (3-1) 

with initial conditions 



0(x,O) = /(x), 
^0(x,O) = fc(x), (3-2) 



we get the formal solution of it. 



r-r „ / X sin ^/At , , , , , 

{x,t) = cos VAtf{x) + ^^k{x). (3-3) 



However in this expression, we need the spectral decomposition of A to define v A explicitly. 
This is not so difficult but the above expression is not useful to obtain the asymptotic 
behavior of the solution, so we want an another expression of the solution. This is done 
by the Duhamel principle. The Duhamel principle is stated as follows. First we define 
G{t) := 6{t) sin y/At/y/A, where 6{t) is the Heaviside function. Then, taking care of the fact 
that the integration kernel G{t; x,y) of G{t) at the region t > is the same as the solution 
of ( |3-1|) with the initial condition 0) = 0, dt4>{x, 0) = 6{x — y), we see that the operator 
G(t) is the solution of the equation, 

—Git) + AG{t) = 5it)l. (3-4) 

From this expression (|3-4| ), the form of G{t) is given by Fourier-transforming G(t) with 
respect to t, 

G{u) = 



-{cu + iOy + A' 
G{t) = J ^Gi^)e-'^'- (3-5) 

Here u+iO means that when we Fourier-transform G{uj) back to G{t), we take the integration 
path so as to deform it into the upper half plane at the spectrum of A. This is because the 
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operator G{t) is only supported in the region t > 0. In Ref. a similar propagator as a 
solution of ( p-4| ) is given, but it is the Feynman propagator obtained by Wick rotation from 
Euclidean space and does not suit our problem. In other words, it obeys different boundary 
condition. So, we get the operator G{t) = 9{t) sm.\/~At/ \/~A^ then we also get the operator 
9{t) co^\/~At by differentiating G{t) with respect to t > 0. Lastly we go back to (pj-3D and 



get the solution of ( |3-1| ) with locally-distributed initial condition. From the fact that the 
kernel G(t; x, y) itself is almost the same as the solution of ( |3-1|) with locally-distributed 
initial conditions (pi^x, 0) = 0, dt4>{x, 0) = 6{x — y), we realize that G(t] x, y) includes all the 
information of asymptotic behavior from the localized disturbance. So below we analyze this 
kernel G{t). In addition to this fact, we also notice that when the localized source term is 
introduced, then from (|3-4| ) we can see its contribution to the wave by simply multiplying 
the kernel G{t) on the source. The Duhamel principle works powerful independently how we 
create a wave. 

Now let us return to our subject to examine the asymptotic behavior of the gravitational 
wave. We formulate our RS2 model so as to use the Duhamel principle For such 

a purpose, we transform hfj_y to and Iz = z,ld = a. Then the equation of motion (|2-5|) 
with boundary condition changes into 



r]^''d,d,<P + {dl - — )0 = 0, 



15 
4^2' 
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+ —(j) = at z = a. (3-6) 

This change of variable hy^y ^ is necessary for achieving the self-adjointness oi A = 
—S'^Widj — d'l + l'b/Az^ in L^-space. Here i,j run from one to three. We choose L^-space as the 
definition domain of A and this matches our physical situation: The AdS space is conformally 
flat and the speed of the signal is finite. Therefore, when we create the disturbance around 
the brane, such a signal does not reach at the boundary of the coordinate system ( p-4|) in 
finite time. So in the finite time observation, the signal is included in the space. We also 
ignore the spin dependence of hpx. According to the Duhamel principle, we get the following 
retarded Green kernel G{uj) by the use of the eigenfunction expansion for this differential 
equation system (p^): 

G{uj]x,z]y,w) 



(27r)3Z(a) _(t^ + ,o)2 + |fc|2 J (27r)3_(^ + ,o)2+ |A?|2 + a' 

Here (t)Q{z) and 0(A, 2;) are the fifth dimensional eigenfunctions and correspond to mass- 
less graviton and the Kaluza-Klein mode respectively, ct(A) is the fifth dimensional spectral 
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measure with support [0, oo) depending on the normahzation of 0(A, z), and Z{a) is a nor- 
mahzation factor of (poiz): 



0(A, 2) = I VA^ (Yi{y/Xa)J2{VXz) - Ji(v^a)r2(v^;2)) , (3-8) 

and Jfj^{z),Y^{z) are the Bessel functions. The derivation of this eigenfunction expansion 
( ]3-7| ) is given in Appendix ^. 

So we can get the analytic expression of the solution of Cauchy problem, then we estimate 
the asymptotic behavior t — oo of this expression on the brane z = w = a. The first term 



of (|3-71 ) has a well known analytic formula so here we estimate the second term Gkk- At the 
brane z = a, the eigenf unctions (f){z; A) become a constant independent of A (really these are 
1) so we can concentrate on the spectral measure cr(A) as calculated below. In the following, 
we can set y = so that the gravitational wave is created at the origin on the brane for 
convenience. 

Our main technique is to use the Paley- Winner-Schwartz theorem which claims that the 
Fourier transformation of a compactly supported infinitely differentiable function decreases 
faster than any power, so that the only contribution which is responsible for the asymptotic 
behavior comes around the singular support^ of the original function. Back to the second 



term of (|3-7| ), and we want to know the asymptotic behavior of the Fourier transformation 
of it, so we must look for the singular support with respect to u. This is easily found, and 
= is the only candidate point which belongs to the singular support of G'kk(^)- Other 
point is holomorphic with respect to u. We also notice that when we approach to u = 
along, for an example, the pure imaginary axis, then the J rf^fcdA-integration around the 
origin := A;^ + A = is convergent by power counting so that the point = is not 
the pole but the branch point of a function G'kk(i^)- So we analyze what branch appear 
after the J (i^fcrfA-integration. This is done as follows. First we change the integrand f{p) 
of the second term of ( |3-7| ) (here we suppress the integration variables except for p) into 
an analytic function F{p) on C \ R+ such that discontinuity over the positive real axis is 
/(p) = F{p + iO) —F{p — iO), and also change the p-integration path into the clockwise curve 
around the real axis. Then we can further deform the p-integration path to Co + Ci + C2 of 
figure |I]. For the integration along the path Cq, uj ^ limit is no longer singular, and the 

*•* A singular support is defined as the closure of the points at which the function is not infinitely 
differentiable. 
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Fig. 1. The contribution to the singular support from the spectral of A. Here uj really means 
Lo + iO which is a pole of F{p). 

singularity of the original f d^fedA-integration concentrates on the p-integration along the 
small closed curve Ci and C2. This last integration is easily performed and if we write the 
integrand in the form F{p) = S{p)/ {~{uj + iO)"^ + p^), then the contribution from the paths 
Ci and C2 is 

- (rH - . (3-9) 

UJ 

Because we search singular support, the holomorphic part of (|3-9| ) with respect to u can 
be ignored. In this way, after performing the integration of remaining variables, we get the 
singularity of Gkk(<-^) as 

4 

GkkH ^«3^aVlogu; + ••• , (3-10) 

OiTT 

in the small u limit. The actual calculus is performed in Appendix |^. Here u really means 
u + iO. We then come back to the time representation of Gkk(^) and get the asymptotic 
behavior 

GKK(t, 0, a; 0, a) ^ --^^(t) + • • • • (3-11) 

If we take account of the Lorentz covariance on the brane, we finally find the Kaluza-Klein 
contribution to asymptotic behavior 

4 

GKKit) 7-5 1^12X3 ' (3-12) 

in the large t limit. Here the suffix + means that the support is only on the t > |x| by 
the causality. This causality could be checked by the \x\ dependent term in ( p-lOp which is 
higher order in u. 

In summary, including also the contribution from the massless graviton mode, the leading 
asymptotic behavior from localized disturbance is 

S(t-\x\) A , , 

2Tia\x\ vr [t^ — 
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The first term is an usual four-dimensional retarded Green function which comes from the 
graviton and represents the Huygens principle in 4-dimensional space-time, that is to say the 
wave form is determined only by its past light-cone. The second term is our main result which 
comes from the continuous spectrum of Kaluza-Klein modes and represents the breaking of 
Huygens principle, that is to say the wave form depends on all of the causal past region. 
This breaking of Huygens principle which arose from the existence of extra dimension, can 
be compared with the flat 5-dimensional space-time's one without brane. In the usual flat 
5-dimension, the retarded Green function behaves as t~^. The difference from our case 
shows the property that in RS2 model, the Kaluza-Klein modes highly decouple from the 
brane mode as noted in Ref. |^). But when we compare this effect with the compact brane 
world case, then in compact case the contribution from the discrete Kaluza-Klein modes 
lead to asymptotically exponential decay. So the power decay behavior of gravitational wave 
distinguishes the noncompact brane world from compact one's. The difference between them 
is caused by whether or not there exists continuous spectrum as seen by our calculus. This 
continuous spectrum, so the power decay of the wave, is a model independent feature of non- 
compact brane world which generically holds even if the continuous spectrum begins after 
some gap (in such a case the simple power decay is replaced by power decay with oscillation.). 
Also this power decay feature is not restricted to the gravitational wave, but extended to any 
other localized mode which comes from the field in the non-compact bulk. So even at the 
classical level and low energy, we can distinguish non-compact brane world from compact 
one (if the power decay behavior different from pure 4- dimensional 's one is observed). Also 
we can conclude that presently observed fields classically like electro-magnetic field, must be 
a non-bulk field living on the brane or must propagate on small compactified space as far 
as they satisfy the Huygens principle in the category of power decay behavior. In the old 
Kaluza-Klein cases without brane, local disturbance in such circumstance cannot be created 
so that the power decay behavior in our calculus distinguishes itself from them. 

So far, we have considered the asymptotic behavior of the gravitational wave which is 
initially distributed around the brane, but then it is interesting to ask whether or not such 
a signal can be observed in a real experiment. As is claimed before, more realistic signal 
depends on how it is created by, for an example, a binary star. We do not pursue such a 
thing, but we estimate the modification of Newton law by the same mode as the above 
contribution. In the process of deriving the asymptotic behavior, the essential thing is that 



the equation ( |3-7| ) with z = w = ais like the Kallen-Lehmann spectral representation except 
for the prescription of the integration path, 

G{uj,x,a;y,a) = / jt^^—, — --ttt^-t^ + / TTrr^dn- 



(27r)3 -(cj + z0)2 + F J (27r)3 ^ -(cu + iO^ + k^ + ^i^'^ 
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and the spectral function behaves in the small /i limit as 

^(^2) ^ ^ ^(27 - 1) + aV log ^ + • • • (3-15) 



Our calculus is concerned with the third term of ( p-15| ), the first and second term does not 



contribute to asymptotic behavior of gravitational wave as seen in Appendix |B[ When we 
compare with the Newton law, the first term mainly contributes and leads to potential 
as is calculated in Ref. ^). Our third term has tinier effect log r/r^ to the gravitational 
potential, so unfortunately it is difficult to observe the signal of RS2 model by the present 
gravitational experiment. 

Finally, we comment on the complementarity between AdS/CFT correspondenceEi)S'0) (see 
a review |l^) for an example) and RS2 geometry. Its original motivation!) is the above 
potential, which is equivalent to the graviton-conformal mode-graviton diagram. Here con- 
formal mode means two point function of stress tensor of dual conformal field theory. This 
interpretation is explicit in the Feynman prescription version of ( |3-14|) and its expansion 



(|3-15|) . As far as we are concerned with the four-momentum dependence, the first term of 



(|3-15|) corresponds to the above mentioned diagram, and the second term represents graviton 



to conformal mode diagram and the third term which is responsible for the asymptotic be- 
havior of gravitational wave, represents graviton to conform mode to graviton to conformal 
mode to graviton diagram (this is k^logk expansion.). So the complementarity between 
AdS / CFT correspondence and RS2 geometry does work qualitatively in our calculus. 

§4. Comments and Discussions 

In this paper, we have considered the Cauchy problem in the RS2 brane world and ana- 
lyzed the asymptotic behavior of the gravitaional wave created by the localized source. The 
power decay behavior resulting from the continuous Kaluza-Klein mode, which is universal 
in noncompact brane world represents the difference (although being difficult to observe it 
now) even at the classical level from a compact brane world, in which case the gravitational 
wave decays exponentially by a mass gap. It also gives a distinction from a Kaluza-Klein 
scenario, in which case we are uniformly distributed over extra dimensions so that no local 
disturbance cannot be created. We also find that in a real calculus, we only need the behavior 
of spectral function of Kaluza-Klein modes around zero or a mass gap. This behavior leads 
to the correction of Newton law or the modification to asymptotic behavior of gravitational 
wave. Especially, in RS2 model, we can get the meaning of spectral function from AdS/CFT 
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correspondence. However at this point, we have a natural question. When we are confined 
in the brane, how can we get the information hke curvature or the shape of extra dimension? 
Such an information is the most important one so that we determine the spectral function. 
Against this question, we have a partial answer in one extra-dimensional case. If the spec- 
tral function is determined, then we get the potential (in our case, the curvature of extra 
dimensions) by solving the Gel'fand-Levitan equation.H' So the remaining question is how 
we observe the full spectral function not restricted to the information around zero. This is a 
typical inverse problem and if we discover the extra dimension in future, it will be important 
to determine the curvature of extra dimension from the signal. Of course, from the point 
of view of inverse problem, more practical method is desired. For example, high frequency 
expansion (eikonal expansion) of gravitational wave uj ^aj{t,x,z) over 

any five-dimensional metric ds"^ = e^^^^^dx^ + dz^) tells us that the leakage rate per unit 



proper time is proportional to m^JGzzi here m means that the Kaluza-Klein mass which 
leaks to fifth dimension and Gzz means the fifth Einstein tensor around the brane. Among 
such things or combination of them, what methods are more appropriate to explore the extra 
dimension? And we ask "Can one feel a shape of extra dimension?" . 
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Appendix A 



Eigenfunction expansion of the operator in ^3- (\ ) 



In this appendix, we will explain the (generalized) eigenfunction expansion of the z-part 
of the differential operator in (|3-6| ) and its completeness . Generalizing the operator a little 
bit, we here consider the self-adjoint differential operator. 



with the boundary condition. 



,p\a) = ^ — -ij{a). (A-2) 



We work in ([a, oo)) space with boundary condition ( [A -21 ), and Eq. ( |3-6| ) corresponds to 
the V = 2 case. Our treatment is standard (see |2T|) for an example) and the outline is as 
follows. First we construct the Green function G{x,y; A) which is the integration kernel of 
the resolvent G{P; A) = (A — P{x, D))^^. Here A e C is included in the complement of the 
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spectrum of P{x,D). This Green function is given by the solutions of differential equation 
P{x, D)tp{x) = Xip{x) as explicitly calculated below. Then using the delta function formula, 



5{x) = — lim ■ 



e^o 27ri \x + ie x — ie^ 
we get the desired form of eigenfunction expansion (the spectral decomposition), 

fix)= [ dX5iX-Pix,D))fix) 
Jn 

= - [ rfA lim — [G(P; X + ie) - G{P; X - ie)] f{x). 

Below we will write this fact explicitly in the form 

S(x — y) = lim / dX \G{x, y; X + ie) — G{x, y; X — ie)] 

2m e^o 

In order to construct the Green function, we consider the differential equation, 

P{x,D)ij{x) = Xip{x), XeC, 
Take two independent solutions of ( |A-6| ) with following initial conditions, 

ipi{a) = cos a, i^'iif^) = sin a, 

■^2(0.) = — sin a, 'V^2('^) = cos a. 

Here a is arbitrary in principle and in the case of ( |A-2| ), we take, 

1/2 - u 



(A-3) 



(A-4) 



(A-5) 



(A-6) 



cos a 



sma 



(A-7) 



(A-8) 



v/a2 + (l/2-z/)2' + (1/2 - u)^ 

Then ipi is the solution of ( |A-6| ) with boundary condition ( |A-2| ). In our case the solutions 
ipi and are given explicitly by Bessel functions, 

-^lix; X) = ^/x AxJuiVXx) + BxY^{VXx) , 

ij2{x; A) = ICxMv^x) + DxY,{v^x)] . (A-9) 
Here A\,B\,Cx, Dx are constants. 



Ax = 


2 


Y,{y/Xa 


Bx = 


1 

2 


Ji,{VXa 


Cx = 


1 

7ra2 
2 


YuiVXa 


Dx = 


1 

2 





1,1 

1(1 _ 

a 2 



cos a (- — u) sma 



\ a 



cos a J — \/XJ,^^i{\^a) cos a 



— VXY,^_i{VXa) si 



sm a 



1.1 



H — ( z/) sina ) + V^Jj^_i(V^a) sma 

a 2 



(A-IO) 
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We also need the fast- decreasing function at infinity for A G C \ R , because we consider 
the L?' space. We write such a function il)'^^'^{x] A) for ImA > 0, and 11)^"^^%] A) for IniA < 0. 
They are given by 



iIj^'^\x; A) = J^(v^x) + iY^{^/Xx 
tlj^'^\x- A) = \jy{VXx) - iY^{V\x 
Using these, we can construct the Green function G^{x, y; A), 
1 



(A-ll) 



G {x,y; A) 



ijji{x)ip^^\y) X < y 



Cx + tPx 
' Ax + iBx 



i^i{x)i)i{y) + i)i{x)i)2{y) 



X < y 
X > y 



for ImA > 0, 



G {x,y]X) 



'ipi{x)^^'^\y) X < y 
iy[^i,^(2)] I ^(2)(x)^i(y) x>y 

Cx - iDx 



^Bx 
iDx 



i'i{.x)ilJi{y) + V^i(a;)V'2(y) 
i'i{.x)ilJi{y) + i'2{.x)'ilJi{y) 



X < y 
X > y 



for ImA < 0. (A- 12) 



Ax - iBx 

Here W[%Ij, 0] = ip{x)(l)'{x) —ip'{x)(f){x) is a Wronskian, so it is a constant and explicitly gives 
Wl'ipi^ip'^^^ = 2i{Ax + iBx)/n and ly V^^^)] ^ ~2i{Ax - iBx)/n. It is easy to check that 
this Green function G^{x, y; A) certainly maps the space with boundary condition (|A-2|) 
to itself by the fact that ipi obeys the boundary condition and ip^^\ip^'^^ are fast-decreasing 
functions for ±ImA > respectively. 

So we get the Green function, then consider the eigenfunction expansion. In such a 
construction, an essential point is that the solution of (|A-6| ) with A-independent boundary 
conditions like (|A-7|) has the holomorphic dependence on A G C. This statement is easily 
checked in our case and also proved in more general setting. And so, to construct the 
eigenfunction expansion using the formula (|A-4|), it is convenient to use the solution ipi,4'2, 



not ip^^\ that is the second expression of in ( |A-12| ). Then we can easily see the 
pole (which corresponds to the discrete spectrum) and the discontinuity across the real axis 
(which corresponds to the continuous spectrum) of the function G^ — G^ with respect to 
A G R because such an information concentrates on the coefficient of ?/'i(a;)^/'i(y) of ( A-12 ). 
Finally we get the delta function in the case of — sin a + 

1 



6{x - y) 



dX 



Al + Bl 



, „ — z/) cos a 7^ 0, 
i)i{x] \)ipi{y] A). 



(A-13) 
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In the case of — sin a + -(o " ^) cos a = 0, we get an extra pole term and this leads to 



oo 



S[x - y) = ^-ilJo{x)ipQ{y) + ^ I dX ^ ipiix; X)ilJi{y; X). (A-14) 

Here Z{a) is a normalization factor and ^0(2^) is a A — limit of ipi{x; A). Their explicit 
forms are given by 



a \ a'^ I J 



—j cos a. (A-15) 

This is the desired form in ( p-7|) in the text. The functions 0o(-2) and 0(A, z) used there are 
iIjq{z)/ cos a and 'ipi{z] A)/ cos a respectively in the case of z/ = 2. 

Appendix B 



The calculus of singularity of Green kernel (\3- Tj ) 



In this appendix, we will give the derivation of (p-10|) which is an expression of singularity 



of ( p-7D . First we change the integration variables in (|3-71 ) to 

^/X = pcos9, k = psm6il2- 

The range of variables are p G [0, 00) , i72 £ and 9 G [0, 7r/2]. Then the the second term 
of ( p-7| ) becomes 

a{^yp cos e)e'P''''^^^-^ cos 9 sin^ 9dpd9dn2 

-(cu + 20)2 + p2 • ^ ■ > 

Here dQ2 is the volume form of 2-sphere. We want to express the integrand of ( |B-2| ) as 
the discontinuity of F{p) = S{p)/ (— (u; + zO)^ + p^) with respect to p. However the full 
form of F{p) and S{p) is not required: As explained in the text, the singularity of ( [B-2| ) is 
concentrated on the expression (p-9| ) and we also know that the higher power of uj is Fourier- 
tansformed to less power of t, so to get asymptotic behavior with respect to t (it is original 
motivation in the text), it is enough to examine the first few singular terms in (|3-9| ), in other 
words the first few terms of ^(p). By definition, the first few terms of ^(p) come from the 
ones of the numerator of integrand in (P-2D . Really the numerator behaves in the small p 
limit as 

a 0^(27-1) 2 2/1 2 2/)i apcos^ 
- + p cos ^ + y P cos 9 log + ■ ■ 



ipsine^?2-af P^QS^sin 9d9dQ2 .^ ^^ 

47^3 
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This expansion corresponds to the one of (p-15|) . From (p-3| ) we easily read the form of S{p). 
However we only need the singularity of ( |3-9| ) so again the full expression is not required, 
for the even function part of ( |B-3| ) does not contribute to the singularity of ( |3-9| ). After all, 
the following expression which discontinuity gives the third term in parenthesis of ( |B-3[ ), is 
needed instead of S{p): 

coi? e {log pf + ^p^ coi? e log p. (B-4) 



Here we take the cut of logarithm along the positive real axis and we ignore the function 
which discontinuity gives the even function part of third term in (P-3|) by the same reason 
above. We make comments on above abandonment of the even function part. In the flat 
4-dimensional case, the same procedure as above only gives the even function part, so no 
power decay behavior exists and it is consistent with the Huygens principle. In the flat 
five dimensional case, the odd function appears so that there exists power decay behavior 
as noted down in the text. Then we go back to the singularity of ( |3-9| ). It is given by 
substituting ( |B-4|) into ( p-9|) with measure in ( |B-3|) : 

• 3 

uo^ logu e^^ '''' cos^ e sin^ 6dddn2 



1677 
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This gives ( p-10 ) in the small uj limit as desired. 



- • 

log u ' " ^ cos^ e sin" OdO. (B-5) 



sin(ci; 


X 


sin 0) 


u 


X 


sin 9 
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